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. Abstract. We study the ionization equilibrium of Fe using photoion- 

I ization models that incorporate improved values for the ionization and 

T-H ■ recombination cross-sections and the charge-exchange rates for the Fe 

I ions. The previously available photoionization models predict concentra- 

^2 ■ tions of Fe^"*" which are a factor of 3-8 higher than the values inferred 

I from emission lines of [Fe ill] and [Fe iv]. Our new models reduce these 

^ ■ discrepancies to factors of 2-5. We discuss the possible reasons behind 

T-H I the remaining discrepancies and present an updated ionization correction 

■ factor for obtaining the Fe abundance from the Fe abundance. 

^: 
CI,: 

2 ■ !• Introduction 

, Iron is expected to be in two main ionization states in H ii regions: Fe"''^ and 

^ ' Fe^^ (with the Fe"*" ion having a non negligible concentration only for very low 

• i-H . excitation nebulae). Therefore, the measurement of Fe iii and Fe iv emission 

^ ' lines will allow us to determine the Fe abundance in these nebulae. [Fe iii] lines, 

although weak, have already been observed in several H ii regions and starburst 
galaxies (e.g. Izotov &: Thuan 1999; Rodriguez 2002), but [Fe iv] lines are much 
weaker and difficult to observe. Therefore, the Fe abundance in H ii regions is 
usually obtained from [Fe iii] lines and an ionization correction factor (ICF), 
obtained from ionization models, to account for the contribution of Fe^"*". For 
example, the relation: 

O ^ 0+ x(Fe++) 0+ ' ^ ' 

where x{X^~^) stands for the ionization fraction of the X"^ ion, is especially 
well suited for determining the Fe abundance from optical observations of H ii 
regions. The ionization potentials of the O and Fe ions are similar (30.6 and 
54.8 eV for Fe++ and Fe^+ , 35.3 and 54.9 eV for 0+ and 0++). Furthermore, 
since both and O^^ can be measured from strong optical lines, one can get 
the O abundance 0/H = 0+/H+ + 0++/H+, and hence also Fe/H from Fe/0. 
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However, available measurements of some weak [Fe iv] lines (Rubin et al. 
1997; Rodriguez 2003) imply Fe^+ abundances which are smaller than the val- 
ues implied by relation (1) by factors 3-8. This discrepancy translates into an 
uncertainty of up to a factor of 6 in the Fe abundances derived for a wide range 
of objects, from the nearby Orion nebula to the low metallicity blue compact 
galaxy SBS 0335—052. Thus, resolving this problem has important implications 
for our understanding of both the evolution of dust in H II regions and the chem- 
ical history of low metallicity dwarf galaxies (see the contribution by Rodriguez 
& Estcban in this volume). 

In order to check whether the discrepancy is due to errors in the ICFs 
predicted by models, we are studying the ionization equilibrium of Fc using 
photoionization models that incorporate recently improved values for all the 
atomic data relevant to the problem. 

2. Results and discussion 

We use the photoionization code NEBULA (Rubin et al. 1991a,b) with the 
following updates: new photoionization cross sections for Fe''', Fe+"*", and 0+ 
(Nahar & Pradhan 1994; Nahar 1996a; Kjeldsen et al. 2002; Verncr et al. 1996), 
new recombination coefficients for Fe"*"^, Fe^"*", O"*" and O"*""*" (Nahar 1996b, 
1997, 1999), all the charge-exchange reactions involving these ions (Kingdon & 
Ferland 1996; and the ORNL/UGA Charge Tranfer Database for Astrophysics 
— http://www-cfade.phy.ornl.gov/astro/ps/data); and the NLTE model stellar 
atmospheres of Pauldrach et al. (2001) with solar metallicity. 

The ionization cross sections are not smoothly varying; they show several 
resonances at different energies. Since the positions of these resonances have 
uncertainties of a few percent, we smoothed (or used available smooth fits) both 
the ionization cross sections and the model atmospheres by convolving with 
Gaussians of widths 3% and 1.5% (in energy), respectively. 

Figure 1 shows the ICFs obtained from different models with T^g = 35000- 
50000 K, total nuclcon density N = 100-10000 cm~^, and Orion metallicity 
Z, Z/10, and Z/30. The results of previous ionization models are also shown 
for comparison. The model results are compared with the results obtained by 
Rodriguez (2003) for several objects with available measurements of [Fe iv] lines. 
The new model results, although still higher than the measured values, are much 
closer to them. The remaining discrepancies could be due to: (1) the need for 
further improvements in the photoionization models, (2) errors in the atomic 
data (in particular the collision strengths) used to derive the Fe^"*" and Fe^^ 
abundances. The effect on the calculated values of a;(0^)/a;(Fe"'"^) of either 
decreasing the Fe"*""*" abundance or increasing the Fe^+ abundance by a factor 
of 2 is also shown in Figure 1. This change would make the results for IC4846, 
30 Doradus (its upper limit), and M42 compatible with the model results. The 
possibility of such a change is provided by the recent calculations of collision 
strengths for Fc++ by McLaughlin et al. (2002). However, they only calculate 
the collision strengths for transitions between terms. Hence, the new collision 
strengths cannot be used either to derive abundances or to check their reliability 
by comparing their predictions with the observed relative intensities of various 
[Fe III] lines. 
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Figure 1. ICFs obtained from our new models compared with the 
previous model results of Stasihska (1990, S90) and Gruenwald & Vie- 
gas (1992, GV92). The squares with error bars show the results ob- 
tained in Rodriguez (2003) for objects with available measurements 
of [Fe iv] lines: from left to right, the planetary nebula IC4846, 
SMC N88A (one position in the object), the low metallicity blue com- 
pact galaxy SBS 0335-052, SMC N88A (a second position), 30 Doradus 
(an upper limit), and M42 (the Orion nebula). The asterisks show the 
effect on the calculated values of either decreasing the Fe^"*" abundance 
or increasing the Fe^+ abundance by a factor of 2. 



Even if the factor of 2 decrease in the Fe"*""*" abundance proved to be real, it 
would not explain the discrepancies shown by the objects with lower metallicities, 
SMC N88A and SBS 0335-052. However, the models with lower metalhcity 
show slightly lower values for the ICF; thus further refinements (such as the 
use of low metallicity model atmospheres) might improve the fit. On the other 
hand, since different [Fe iv] lines were used to derive the Fe^+ abundance for 
the five objects, changes in the relevant collision strengths might explain all the 
discrepancies. 

In any case, the Fe abundances implied by [Fe III] emission lines should be 
obtained through the relation implied by a fit through the new model results in 
Figure 1: 

Fe _ f 0+ y-'' Fe++ 

o " 0++ 0+ • 
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